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ABSTRACT 

We present an analysis of UBVRcIcJH photometry and phase-resolved optical spec- 
troscopy of NSVS 14256825, an HW Vir type binary. The members of this class consist 
of a hot subdwarf and a main-sequence low-mass star in a close orbit (Porb ^ 0.1 d). 
Using the primary-eclipse timings, we refine the ephemeris for the system, which has 
an orbital period of 0.11037 d. From the spectroscopic data analysis, we derive the 
effective temperature, Ti — 40000 ± 500 K, the surface gravity, log 51 = 5.50 ± 0.05, 
and the helium abundance, n(He)/n(H) = 0.003 ± 0.001, for the hot component. 
Simultaneously modelling the photometric and spectroscopic data using the Wilson- 
Devinney code, we obtain the geometrical and physical parameters of NSVS 14256825. 
Using the fitted orbital inclination and mass ratio {i = 82? 5 ± 0?3 and q = M2/M1 = 
0.260±G.012, respectively), the components of the system have Mi = 0.419±0.070 M©, 
i?i = 0.188 ± 0.010 Rq, M2 = 0.109 ± 0.023 Mq, and i?2 = 0.162 ± 0.008 Rq. From 
its spectral characteristics, the hot star is classified as an sdOB star. 

Key words: Binaries: eclipsing - stars: fundamental parameters - stars: subdwarf - 
stars: low-mass - stars: individual: NSVS 14256825. 
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1 INTRODUCTION 

HW Virginis (HW Vir) systems consist of a subdwarf B or 
OB (sdB or sdOB; hereafter referred to as sdB) plus a main 
sequence star, which form an eclipsing pair in a compact or- 
bit, Porb ~ 0.1 d (Heber 2009 1. These systems are believed to 



evolve through a common envelope phase when the primary 
(sdB) is a red giant. During this stage, the secondary star 
(dM) spirals in towards the primary one and the potential 
gravitational energy released is absorbed by the envelope. 



which is subsequently ejected (Taam & Ricker 20101. The 



final separation between the dM and the sdB stars depends 
on the initial mass ratio of the binary and the initial sepa- 
ration. 

sdB stars consist of a helium-burning core covered by a 
thin hydrogen-dominated envelope. The atmosphere abun- 
dance is normally n(He)/n(H) ~ 0.01, the effective tempera- 
tures are in the range of 22000-37000 K, and the logarithmic 
surface gravities are normally between 5.2 and 5.7. They 
populate a narrow strip on the extreme horizontal branch 
(EHB) in the Hertzsprung-Russell diagram. The single-star 
stellar evolution predicts a narrow mass range; 0.46 - 0.50 



* Based on observations carried out at the Observatorio do Pico 
dos Dias (OPD/LNA), in Brazil 
f E-mail: leonardo@das.inpe.br 



M0 (Dorman et al. 1993). On the other hand, models based 
on binary evolution predict a broader range of masses, from 
0.3 to 0.8 Mq (Han et al. 2003). Hence an important step in 
understanding the origin of sdB stars is the determination 
of their mass distribution. A recent review of sdB stars is 
presented by |Heber| ( |2009[ ) . 



There are currently ten members of the HW Vir class, 
whose main features are summarised in Table IT] Among 
them, NSVS 14256825 (2MASS J2020-H0437; hereafter re- 
ferred to as NSVS 1425) is one of the least studied. It was 
discovered in the public data from the Northern Sky Vari- 
ability Survey ( Wozniak et al.|[2004 |. The sole information 
on this system comes from the photometric data by |Wils,| 
Giorgio & Sebastian (20071. These authors obtained B, V, 



and Ic photometric light curves. The main parameters ob- 
tained by |Wils et al] ( |2007[ ) are listed in Table [l] 



In this paper, we report on multi-band photometry 
and phase-resolved optical spectroscopy of NSVS 1425. We 
present an improved solution for its geometrical and phys- 
ical parameters and discuss these results in the context of 
HW Vir systems and their evolution. 
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Table 1. The currently known sdB+dM eclipsing binaries. 



Name 



Ti 

(K) 



Ml 
(Mq) 



M2 

(Rq) 



Ri 

(Rq) 



R2 

(Rq) 



log gi Period 

(d) 



Refs. 



Notes 



AA Dor 

NSVS 14256825 

PG 1336-018 

2M 1533+3759 
2M 1938+4603 
HS 0705+6700 
PG 1241-084 
HS 2231+2441 
SDSSJ0820+0008 
BUL SC16 335 



42000 

40000 

31300 
32740 
30400 
29564 
28800 
28488 
28370 
26700 



0.33/0.47 

0.25 

0.419 

0.46 

0.389 

0.459 

0.377 

0.48 

0.48 

0.48 

0.47/0.499 

0.25/0.47 



0.064/0.079 

0.054 

0.109 

0.21 

0.110 

0.113 

0.12 

0.13 

0.14 

0.075/0.072 

0.068/0.045 



0.179/0.20 

0.165 

0.190 

0.150 

0.166 

0.230 
0.176 
0.250 



0.097/0.108 

0.089 

0.151 

0.160 

0.152 

0.186 
0.180 
0.127 



5.46 


0.261 


1,2,3 


15,16 






4 


16 


5.50 


0.11037 


this work 


15,16 


- 


0.11040 


5 


15 


5.60 


0.10102 


6 


15,16 


_ 


_ 


7 


17 


5.58 


0.16177 


8 


15,16 


5.425 


0.12576 


9 


15,16 


5.40 


0.09565 


10 


15,16 


5.63 


0.11676 


11 


15 


5.39 


0.11059 


12 


15,16 


5.48 


0.096 


13 


15,16 


- 


0.12505 


14 


15 



Heberl | |2009|l; ^Klepp fc 



2009 



1999 



'IWils 



RauchI | |201l|l; ^Hilditch et al.| ||2003|l; JRucinskil ( |2009 1 
^ jCharpinet et al.| j2008l l ; ^fbr et al.| (|2010t; ^0stensen et al.| | |2010[ l; '^Drechsel 
; l ^0stensenet al.| ( |2007^ ^'^°i°'-°t al Monnli- J^' ' "" 



[Wood fc Saffer 
^Geier et al.|(|2011^; ^^Polubek et al.|||2007|l; ^'^Light Curves; ^''Spectroscopy; ^''Asteroseismology 



et ah| ||2007| l; 
et al.| |200i7 ; 



' Vuckovic et al. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Optical and near infrared photometry 

The observations were carried out using the Observatorio 
do Pico dos Dias (OPD/LNA) facilities, in Brazil. Photo- 
metric data in the U, B, V, Re, Ic, J, and H bands were 
obtained from July to November, 2010. Optical observations 
were performed using a CCD camera attached to the 0.6-m 
TAG telescope. Near infrared photometry data were collected 
by means of a CamlV imager attached to the 1.6-m Perkin- 
Elmer telescope. 

The procedure to remove undesired effects from the 
CCD data included obtaining ~100 bias frames and ~30 
dome flat-field images for each night of observations. The 
NIR flat-field images were produced using separate se- 
quences of 30 "on" and 30 "off' exposures. The resulting 
of the "on" image minus "off' image was used as the mas- 
ter flat-field image. Table [2] summarises the characteristics 
of the data collected for NSVS 1425. In this table, iV is the 
number of individual images obtained with the integration 
time texp. 

The preparation of the CCD data was performed us- 
ing standard IRAHj tasks (Tody, 1993) and consisted of sub- 
tracting a master median bias image from each program im- 
age, and then dividing the result by a normalised flat-field 
frame. In the J and H bands, additional steps of lineariza- 
tion and sky subtraction from dithered images were used in 
the preparation of the data. For both optical and infrared 
data, differential photometry was used to obtain the relative 
fluxes between the target and a set of constant flux stars in 
the fleld of view. As the NSVS 1425 field is not crowded, the 
extraction of the fluxes was carried out using aperture pho- 
tometry. Figure [1] shows a finding chart for NSVS 1425 with 
a reference star and a number of additional comparison stars 
adopted for differential photometry. Photometric standard 
stars were observed each night in order to calibrate the op- 
tical data. The NIR calibration is directly provided by the 



http://www.iraf.noao.edu 
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Figure 1. Finding chart for NSVS 14256825 in the Re band ob- 
tained using the OPD/LNA 0.6-in telescope. The circle shows 
NSVS 1425, the diamond is the adopted reference star, and the 
squares outline additional comparison stars. 



2MASS catalogue magnitudes for the reference and compar- 
ison stars. Figure [2] shows a sample of calibrated light curves 
folded on the NSVS 1425 orbital period. 

As can be seen in Figure [2J the light curves of 
NSVS 1425 show a prominent reflection effect, which in- 
creases towards longer wavelengths. The depth of the pri- 
mary eclipse is ~0.7mag and does not change signiflcantly 
with wavelength, while the depth of the secondary eclipse 
increases towards longer wavelengths, from ~0.1 mag in the 
U band to ~0.18mag in the H band. Table|3] shows the ap- 
parent magnitudes for NSVS 1425 at primary and secondary 
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Figure 2. Calibrated NSVS 1425 light curves in the U, B, V, Re, 
Ic, J, and H bands folded on the 0.1104 d orbital period. The 
Rc-band curve was displaced upwards by 0.5 mag to improve 
visualisation. 



Table 2. Log of the photometric observations. 



UT Date 


N 


tcxpiSj 


Telescope 


Filter 


2010 Jul 30 


300 


20 


0.6-m 


Rc 


2010 Jul 30 


185 


30 


1.6-m 


J 


2010 Jul 31 


250 


30 


0.6-m 


B 


2010 Jul 31 


450 


20 


0.6-m 


Rc 


2010 Jul 31 


177 


30 


1.6-m 


J 


2010 Jul 31 


235 


30 


1.6-m 


H 


2010 Aug 01 


235 


40 


0.6-m 


U 


2010 Aug 01 


75 


40 


0.6-m 


B 


2010 Aug 01 


100 


30 


1.6-m 


Y 


2010 Aug 01 


215 


30 


1.6-m 


J 


2010 Aug 02 


230 


40 


0.6-m 


U 


2010 Aug 02 


400 


20 


0.6-m 


Ic 


2010 Aug 06 


120 


30 


0.6-m 


V 


2010 Aug 07 


220 


30 


0.6-m 


V 


2010 Aug 08 


220 


30 


0.6-m 


Ic 


2010 Aug 09 


160 


25 


0.6-m 


Rc 


2010 Aug 10 


80 


25 


0.6-m 


Rc 


2010 Aug 18 


800 


10 


0.6-m 


Rc 


2010 Aug 18 


530 


15 


0.6-m 


Ic 


2010 Aug 19 


225 


40 


0.6-m 


u 


2010 Aug 20 


420 


20 


0.6-m 


B 


2010 Sep 01 


160 


30 


0.6-m 


V 


2010 Sep 03 


250 


30 


1.6-m 


J 


2010 Sep 04 


142 


50 


1.6-m 


J 


2010 Nov 03 


352 


20 


0.6-m 


Ic 




4000 4200 4400 4600 4B00 
A (A) 

Figure 3. Upper panel: a normalised individual spectrum of 
NSVS 1425 with 900 s integration time. Lower panel: the aver- 
age of 36 spectra after correcting for orbital motion. 



2.2 Optical spectroscopy 

The spectroscopic observations were performed using the 
Cassegrain spectrograph attached to the 1.6-m telescope at 
OPD/LNA. Thirty six (36) spectra were obtained using the 
1200 1/mm grating and integration times of 10 or 15 minutes. 
The spectral coverage of this configuration is 3950-4900 A, 
using 1.8 A resolution (from the FWHM of the wavelength 
calibration lines). A hundred (100) bias frames and 30 flat- 
field frames were obtained each night to remove systematic 
signatures from the CCD detector. Observations of a He- 
Ar comparison lamp were made every two exposures on the 
target to provide wavelength calibration. The spectropho- 
tometric standard s tars HR 1544, HR 7596, and HR 9087 
( jHamuy et al.|[T992[ ) were observed for flux calibration. 

The reduction of the spectra was carried out following 
the steps of bias subtraction, flat-field structure removal, 
optimal extraction, wavelength calibration, and flux cali- 
bration using the standard routines in IRAF. The upper 
panel in Figure[3]shows a typical normalised individual spec- 
trum. The lower panel shows the average of all spectra after 
Doppler shifting using the radial velocity orbital solution 



(see Section 3.2 1 



3 ANALYSIS AND RESULTS 

3.1 Ephemeris 

To determine an ephemeris for the times of the primary min- 
imum in NSVS 1425, we combined our timings and those 
from |Wils et al. | ( |2007| l, after converting them to barycen- 
tric dynamical time (TDB). Our eclipse timings were ob- 
tained by modelling the primary eclipse using the Wilson- 
Devinney (WD) code ( [Wilson fc Devinney||197l] > together 
with a Markov chain Monte Carlo (MCMC) procedure 
(Gilks, Richardson & Spiegelhalter||1996l to obtain the un- 



certainties. The geometrical and physical parameters of the 
system are calculated as described in Section [3. 4[ These pa- 



4 L. A. Almeida et al. 



Table 3. Apparent magnitudes of NSVS 1425 at primary and 
secondary minima. 



UT Date 


Band 


Primary 


Secondary 






mmimum 


mmimum 


2010 Aug 02 


U 


12.41±0.20 


11.73±0.20 


2010 Aug 20 


B 


13.76±0.15 


13.06±0.15 


2010 Aug 07 


V 


13.93±0.13 


13.23±0.13 


2010 Jul 31 


Re 


13.79±0.12 


13.10±0.12 


2010 Aug 18 


Ic 


14.18±0.16 


13.47±0.16 


2010 Aug 01 


J 


14.50±0.24 


13.80±0.24 


2010 Jul 31 


H 


14.60±0.25 


13.91±0.25 



rameters were used as fixed inputs for the WD code and 
only the times of individual primary eclipses are left as free 
parameters. The median values and the 1-cr uncertainties ob- 
tained from the marginal distribution of the fitted instants 
of minimum were adopted as the best values for location and 
error of each timing. Using the expression Tmin = Tq + ExP, 
where Tmin are the predicted times of primary minimum, To 
is a fiducial epoch, E is the cycle count from To, and P is the 
binary orbital period, we obtained the following ephemeris, 

Tmin = TDB 2454274.2087(1) + 0.110374230(2) x E. (1) 



3.2 Radial velocity solution 

The radial velocities were obtained using the task FXCOR in 
IRAF. Initially, a combination of all 36 spectra was used as a 
template for the cross correlation with the individual spec- 
tra. Regions around H/3, H7, H(5, He, and Hell A4686 were se- 
lected (see Figure|3| to improve the signal/noise ratio of the 
correlation procedure. The resulting radial velocity solution 
was used to Doppler shift all individual spectra to the or- 
bital rest frame. A better quality template is then produced 
from these rest-frame spectra. The procedure was iterated a 
number of times until the radial velocity solution converged. 
Table HI lists individual radial velocities and Figure |4] shows 
the radial velocity curve folded on the orbital phase together 
with the best solution for a circular orbit. The modelling 
provides a semi-amplitude Ki — 73.4 ± 2.0 km s~^ and a 
systemic velocity 7 = —12.1 ± 1.5 km s^^. 



3.3 Atmospheric parameters 

The atmospheric parameters of the sdB star can be deter- 
mined using the Balmer and helium lines in the blue range 
of the spectrum. The spectra obtained in a 0.1 phase in- 
terval centred in the secondary eclipse, i.e., 0.45-0.55, were 
used to minimise the contribution of the refiection effect. 
Using x^ as the figure of merit, the combined spectrum was 
matched to a grid of synthetic spectra retrieved from the web 
page of TheoSSAr] The synthetic spectra were calculated by 
the Tiibingen non-local thermodynamic equilibrium Model- 
Atmosphere Package]^ (TMAP). Two different metallicities 
were used: Model A with zero metallicity; and Model B with 



the metallicity adopted by Klepp & Ranch (2011 1. The grid 
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Figure 4. Radial velocity of the prominent lines in the spec- 
tra of NSVS 1425. The phases are calculated using the ephemeris 
presented in Section |3.1| 



Table 4. Spectroscopic observations and radial velocities. 



UT Date 


BJD{TDB) 


tcxp 


V 


Orbital 




(2450000+) 


(s) 


(Ijms-l) 


Phase 


2010 Set 01 


5441.48256 


900 


2.29±12.29 


0.60 


2010 Sot 01 


5441.49251 


600 


-8.05±11.86 


0.69 


2010 Sot 01 


5441.50371 


900 


66.46±10.50 


0.79 


2010 Set 01 


5441.51799 


900 


12.46±10.60 


0.92 


2010 Sot 01 


5441.52904 


900 


-40.38±12.10 


0.02 


2010 Set 01 


5441.54078 


900 


-67.43±10.21 


0.13 


2010 Sot 01 


5441.55825 


900 


-83.52±11.16 


0.29 


2010 Set 01 


5441.56970 


900 


-64.55±10.12 


0.39 


2010 Set 01 


5441.58084 


900 


-21.99±12.32 


0.49 


2010 Set 01 


5441.59792 


900 


53.03±12.84 


0.65 


2010 Set 01 


5441.60920 


900 


60.29±10.98 


0.75 


2010 Set 01 


5441.62057 


900 


47.18±9.33 


0.85 


2010 Set 01 


5441.63857 


900 


-44.83±10.33 


0.02 


2010 Set 01 


5441.65259 


900 


-78.09±12.25 


0.14 


2010 Set 01 


5441.66364 


900 


-76.39±12.18 


0.24 


2010 Set 02 


5442.42478 


900 


-63.45±10.09 


0.14 


2010 Set 02 


5442.43581 


900 


-84.04±10.91 


0.24 


2010 Set 02 


5442.45306 


900 


-37.19±9.84 


0.40 


2010 Set 02 


5442.46413 


900 


-6.20±8.10 


0.50 


2010 Set 02 


5442.48116 


900 


44.56±8.03 


0.65 


2010 Set 02 


5442.49223 


900 


65.25±9.62 


0.75 


2010 Set 02 


5442.51304 


900 


32.99±10.88 


0.94 


2010 Set 02 


5442.52404 


900 


-32.57±8.55 


0.04 


2010 Set 02 


5442.54110 


900 


-82.53±7.39 


0.19 


2010 Set 02 


5442.55209 


900 


-77.39±8.10 


0.29 


2010 Set 02 


5442.56870 


900 


-30.80±10.13 


0.44 


2010 Set 02 


5442.57970 


900 


20.41±10.04 


0.54 


2010 Set 02 


5442.60086 


900 


57.22±11.46 


0.74 


2010 Set 02 


5442.61186 


900 


58.08±9.66 


0.83 


2010 Set 02 


5442.62893 


900 


-8.53±10.40 


0.98 


2010 Set 02 


5442.63992 


900 


-65.09±9.83 


0.09 


2010 Set 02 


5442.65664 


900 


-91.21±10.98 


0.24 


2010 Set 02 


5442.66763 


900 


-70.64±11.55 


0.34 


2010 Set 02 


5442.68442 


900 


-3.21±6.66 


0.49 


2010 Set 02 


5442.69542 


900 


23.22±10.57 


0.59 


2010 Set 02 


5442.71246 


900 


47.40±9.12 


0.74 



^ http://dc.g-vo.org/theossa 

3 http://astro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 



is composed by: 26 values of effective temperatures, 30000 
K < T ^ 43000 K with 500 K steps; 16 surface gravities, 
5.2 ^ logg ^ 6.0 with 0.05 dex steps; and 10 helium abun- 
dances, 0.001 ^ n(He)/n(H) ^ 0.01 with 0.001 dex steps. 
All synthetic spectra were convolved with the projected ro- 
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Figure 5. The best fits to the Balmer and heUum fines used to 
derive effective temperature, surface gravity, and hefiuni abun- 
dance. The observed spectra are presented using black dots and 
the red fines represent the best-fit synthetic spectra. The left and 
right panels correspond to models with zero nietallicity and the 
metallicity from |Klepp fc Rauch| | [2011| , respectively. 

tational velocity iisin i = 73.4 km s~^ and with an instru- 
mental profile of FWHM 1.8 A. 

The synthetic Balmer (H/3 to He) and helium 
(HeIA4471 and HeIIA4686) lines were used in the fitting 
procedure to determine the effective temperature, surface 
gravity, and He abundance. The best fit yields T — 38000 ± 
500 K, logg = 5.2 ± 0.05, and n(He)/n(H) = 0.003 ± 0.001 
for the Model A, and T = 40000±500 K, logg = 5.50±0.05, 
and n(He)/n(H) = 0.003±0.001 for Model B. Figure[5]shows 
the observed spectra together with the best synthetic spec- 
tra for both models. The x^ oi model B is ~5 per cent better 
than that of model A. 



3.4 Simultaneous modelling of light curves and 
radial velocities 

In order to obtain the geometrical and physical parameters 
for NSVS 1425, we simultaneously analysed the U, B, V, Re, 
Ic, J, and H light curves and the radial velocity curve us- 
ing the WD code. The original WD code uses a differential 
correction method for improving an initial solution. This 
method works well if the initial parameter values are close 
to those corresponding to the optimal solution. However, if 
they are close to a local minimum, the differential correc- 
tion procedure may fail to find the best solution. To solve 
this problem, the WD code was used as a "function" to be 



optimised by the genetic algorithm PIKAIA ( Charbonneau 



1995), which is adequate to search for a global minimum in 



a model involving a large set of parameters. 

To examine the marginal distribution of probability of 



the parameters and to establish realistic uncertainties, we 
used the solution obtained by PIKAIA as an input to an 
MCMC procedure. 

Due to the large number of parameters to be fitted, it 
is important to constrain them using theoretical and spec- 
troscopic information. From the spectroscopic analysis, we 
adopted the effective temperature of the primary star as an 
initial value. We can also constrain the mass ratio, q, us- 
ing the mass function (Eq. pi and assuming that the mass 



of the sdB star. Mi, 



et al. 1998 Han et al 



semi-amplitude, Ki 



is in the range 0.1 
2003] ), 



0.8 Mq priebel 



and that the radial velocity 
is 73.4 kms"^ 



Ml X ((J sin i)^ 
(I + OP 



= 1.0361 X 10"^(1 - e'^f^'^KfP. 



(2) 



As the components are in close orbit, the timescales of 
synchronisation and circularisation are much shorter than 
the helium burning lifetime ( |Zahn|[T977[ ). Thus, the orbit 
can be considered circular (e = 0) and the rotation of the 
components synchronised with the orbit. Finally, adopting 
the range of orbital inclinations for eclipsing binaries, 75° < 
i < 90°, we obtained 0.21 < q < 0.45 for the mass ratio 
range. 

Mode 2 of the WD code, which sets no constraints on 
the Roche configuration, was used. The luminosity of the 
secondary star was computed assuming stellar atmosphere 
radiation. Linear limb darkening coefficients, Xi, were used 
for both stars. Regarding the sdB star, we used the co- 
efficients c alculated by |Diaz-Cordoves, Claret fc Gimenez] 
( [1995[ ) and [Claret, Diaz-Cordoves fc Gimenez ( 1995 1 for a 
star with effective temperature T = 40000 K and surface 
gravity logg = 5.0. These are the closest values to those of 
the hot component in NSVS 1425 for which limb darkening 
coefficients have been published in literature. On the other 
hand, the linear limb darkening coefficients of the cool star 
were left as free parameters, since the proximity of the hot 
star can significantly change these coefficients with respect 
to those of a single star. As the sdB star has a radiative en- 
velope, its gravity darkening exponent, /3i, and its bolomet- 
ric albedo for refiective heating and re-radiation, Ai, were 
set to unity (Rafert fc Twigg 19801. The gravity darken- 



ing exponent of the secondary component, /32, was fixed at 
0.3, which is appropriate for convective stars (Lucy 19671. 



(20101, Kilkenny et al. 



(19981, and 



As shown by For et al. 
[Drechsel et al.| ( |2001[ ), the albedo of the secondary star, A2, 
can assume physically unrealistic values A2 > 1, especially 
at longer wavelengths where the refiected-reradiated light 
is more intense. For this reason, it was decided to perform 
two modellings: in Model 1, we adopt a constant (but free 
parameter) secondary albedo for all photometric bands; in 
Model 2, we consider variable and independent albedos for 
all photometric bands. 

In both cases the remaining fitted parameters consist of: 
the mass ratio, q = M2/M1; the orbital inclination, i; the 
separation between the components, a; the Roche potentials, 
f7i and 5^2; and the effective temperatures of the two stars, 
Ti and Tj. 

In order to optimise the computational time, all light 
curves were binned with 160 seconds time resolution and the 
error of the bin average outside of the eclipses was assumed 
as the uncertainty. To test the goodness of fit, we use the 
reduced x^ defined as 
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Table 5. System parameters of the best fit for photometric light 
curves in U, B, V, Re, Icj J ^^'^ H bands and radial velocity data 
for NSVS 1425. 



Orbital Phase 



Figure 6. The best simultaneous fit for the U, B, V, Re, Ic, J, 
H light curves and radial velocity data using the WD code. The 
light curves are the same as those shown in Figure [2] displaced 
vertically for better visualisation. 



2 
X = 



1 y^/O, - 



a 



(3) 



where Oi are the observed points, d are the corresponding 
model, di are the uncertainties at each point, and n is the 
number of points. Figure [6] shows the best fit together with 
the hght and radial velocity curves while Table [5] lists the 
fitted and fixed parameters. 

The values found for almost all parameters in the two 
models are consistent. However, the y^ for Model 2 is ~40 
per cent better than for Model 1. 



3.5 Fundamental Parameters 

Physical and geometrical parameters such as masses, radii, 
and separation between the two components of the system 
can be derived from the solutions obtained in the previ- 
ous sections. Substituting the values of the orbital period 
(forb = 0.110374230 d), semi-amplitude of the radial ve- 
locity (K\ = 73.4 km s^^), mass ratio (5 — 0.260), and 
inclination (i = 82?5) in Eq. [2] A/i = 0.419 ± 0.070 Mq 
is obtained for the primary mass. The primary mass and 
the mass ratio are used to derive the secondary mass, 
M2 = 0.109 ± 0.023 Mq. Using Kepler's Third Law, one 



Parameter 


Model 1 


Model 2 


Fitted Parameters 


q = (M2/M1) 


0.28 ±0.013 


0.26 ±0.012 


in 


82.5 ±0.4 


82.5 ±0.3 


n- 


4.58 ±0.13 


4.55 ±0.11 


n- 


2.80 ±0.15 


2.69 ±0.12 


Ti{K) 


42300 ± 400 


42000 ± 500 


T2{K) 


2400 ± 600 


2550 ± 550 


a^ (Rq) 


0.74 ±0.04 


0.80 ± 0.04 


A=(U) 


1.50 ±0.11 


2.0 ±0.15 


^i(B) 


1.50 ±0.11 


1.35 ±0.13 


^^(V) 


1.50 ±0.11 


1.20 ±0.12 


A|(Rc) 


1.50 ±0.11 


1.05 ±0.09 


Alilc) 


1.50 ±0.11 


1.3 ±0.12 


Am 


1.50 ±0.11 


0.95 ±0.14 


^^(H) 


1.50 ±0.11 


1.10 ±0.15 


X2{U) 


0.64 ±0.04 


0.68 ± 0.04 


X2{B) 


0.69 ± 0.04 


0.74 ±0.05 


X2(V) 


0.78 ± 0.03 


0.80 ± 0.04 


X2(Rc) 


0.83 ± 0.02 


0.87 ±0.02 


^2{Ic) 


0.90 ± 0.03 


0.92 ± 0.03 


X2{J) 


0.93 ± 0.04 


0.95 ± 0.04 


X2{H) 


0.98 ± 0.05 


0.99 ± 0.05 


Roche radii'' 


ri (pole) 


0.231±0.006 


0.233±0.005 


Ti (side) 


0.233±0.006 


0.235±0.005 


ri (point) 


0.235±0.007 


0.236±0.006 


ri (back) 


0.234±0.007 


0.236±0.006 


r2 (pole) 


0.180±0.016 


0.194±0.014 


r2 (side) 


0.182±0.016 


0.198±0.016 


r2 (point) 


0.191±0.019 


0.210±0.019 


r2 (back) 


0.189±0.019 


0.207±0.018 


Fixed Parameters 


Al 


1.0 


1.0 


% 


1.0 


1.0 


0.3 


0.3 


xl(\]) 


0.242 


0.242 


x^(B) 


0.233 


0.233 


x^(V) 


0.209 


0.209 


4 (Re) 
4(Ic) 

x((J) 
x((H) 


0.176 


0.176 


0.147 


0.147 


0.112 


0.112 


0.095 


0.095 


Goodness of fit 


X^ 


2.1 


1.2 



" Roche surface potential; 

^ Components separation; 

■^ Bolometric albedo; 

"^ Gravity darkening exponent; 

'^ Linear limb darkening coefficient from |Dia2-Cordoves et al.| 

| 1995| l; 

^ Linear limb darkening coefficient from [Claret et al.| | |l995| ; 
'' In units of orbital separation. 
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Table 6. Fundamental parameters for NSVS 1425. 



Parameter 


Model 1 


Model 2 


Ml (Mq) 


0.346±0.079 


0.419±0.070 


M2 (M©) 


0.097±0.028 


0.109±0.023 


Ri (R©) 


0.173±0.010 


0.188±0.010 


R2 (R©) 


0.137±0.008 


0.162±0.008 


Ti(K) 


42300±500 


42000±400 


T2(K) 


2400±500 


2550±500 


logSi 


5.50±0.14 


5.51±0.11 


log 92 


5.15±0.16 


5.05±0.13 


a{R0) 


0.74±0.04 


0.80±0.04 



can obtain the orbital separation, a = 0.80 ± 0.04 R©, from 
which the absolute radii follow, Ri — 0.188 ± 0.010 R© and 
R2 = 0.162 ±0.008 R©. In Table|6]we show the fundamental 
parameters for NSVS 1425 derived from Model 1 and Model 
2. 



3.6 Rossiter-McLaughlin effect 

An interesting spectroscopic signature in eclipsing binary 



systems is the Rossiter-McLaughlin (RM) effect (Rossiter 
1924[ iMcLaughliiil |1924[). This effect occurs when the 



eclipsed object rotates. There is evidence of this effect on 
the radial velocity curve in the phase interval 0.95 — 1.05 
(Figure [m. Unfortunately, our data are not enough to model 
the RM effect to derive the alignment rotational parameters 
of the NSVS 1425 stars. However, our observed points are 
consistent with the predicted RM effect for aligned rotating 
stars obtained by the WD code with the parameters shown 
in Tabled 
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Figure 7. Position on the (Tcfi, logg) diagram of the hot com- 
ponent of NSVS 1425 compared with other sdB stars (see Ta- 
bloll}. Isolated sdB and sdOB stars presented in |Edel mann (2003|l 
are shown with open circles. Dashed lines represent evolutionary 
tracks for different masses in the post-EHB evolution | |Dorman et| 
|al.|1993^ . The zero age extreme horizontal branch (ZAEHB) and 
terminal age extreme horizontal branch (TAEHB) are represented 
by solid lines. 



4.2 Evolution 



Han et al. (2003 1, from a detailed binary population synthe- 



sis study, presented three possible channels for forming sdB 
stars: 



4 DISCUSSION 

4.1 Characteristics of the primary star 

Of all HW Vir systems, the primary of NSVS 1425 has the 
second highest temperature (see Table [T]) , consistent with 
the prominent HeIIA4686 line (see Figure^. Accordingly, 
we suggest that the primary of NSVS 1425 is an sdBO star 
which means this system is very similar to AA Dor and, 
hence, in a rare evolutionary stage ( Heber|2009[ ). 

Comparing the values of loggr derived from the simul- 
taneous fit to photometric and spectroscopic data (Table |6| 
with those obtained from the modelling of the spectral lines 



of the primary star (Section 3.31, it is clear that the model 



with metallicity equal to that adopted by |Klepp fc Ranch] 
(20111 provides consistent results, whereas the model with 



zero metallicity has a discrepancy. We noticed that the same 
kind of discrepancy had been found by |Rauch| (|2000[ ) in the 
analysis of the primary in AA Dor. Ranch ( 2000[ ) obtained 
logg = 5.21 from spectroscopic data, whereas |Hilditch, Har-| 
ries, fc Hill| ( [T996| ) had derived logy = 5.53 from photomet- 
ric data modelling. This discrepancy was solved by |Klepp fc| 
Rauch (20111 with the improvement of the Stark broaden- 



ing modelling, the minimisation of the reflection effect, and 
adoption of metal-line blanketing. 



(i) One or two phases of common envelope (CE) evolu- 
tion; 

(ii) One or two stable Roche lobe overflows; 
(iii) A merger of two He-core white dwarfs. 



Driebe et al. ( 1998 1 and Heber et al. ( 2003 1 suggest an- 



other scenario to form an sdB star with low mass in binaries 
called post-RGB. This scenario is similar to the channel (i) 
proposed by Han et al. ( 2003 1, except that the resultant sdB 



in the first CE phase has insufficient mass in its core to ig- 
nite helium. However, it will evolve as a helium star through 
the sdB star region to form a helium core white dwarf. 

In FigurelTlwe compare the position of the primary com- 
ponent of NSVS 1425 on the (Teff, log g) diagram with other 
sdB stars in short-period binary systems (see Table [I]). We 
also show a sample of single sdB and sdOB stars analysed by 



Edelmann (20031. In the same graph we display evolution- 



ary tracks for different masses in the post-EHB calculated by 
[Dorman, Rood fc 0'Connell| ( |1993[ ) for single star evolution. 
As can be seen, the position of the NSVS 1425 primary is 
close to that of the AA Dor primary. The evolutionary track 
for single star evolution would only marginally explain the 
mass obtained for the sdB star in NSVS 1425. Out of the 
possible channels to form an sdB in binaries, we suggest 
that the channel (i) is probably the evolutionary scenario 
for NSVS 1425. 
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5 CONCLUSION 

We present a photometric and spectroscopic analysis of 
the NSVS 1425 system. With a short orbital period, P = 
0.110374230(2) d, this binary shows both primary and sec- 
ondary eclipses and a prominent reflection effect. From the 
spectroscopic analysis we obtain 73.4 ± 2.0 km s~^ for the 
semi-amplitude of the radial velocity and — 12.1±1.5 km s~^ 
for the systemic velocity. The atmospheric parameters of the 
primary component (sdOB star), namely, effective tempera- 
ture, T = 40000 ± 500 K, surface gravity, log g = 5.5 ± 0.05, 
and Helium abundance, n(He)/n(H) = 0.003 ± 0.001, were 
calculated matching the observed spectrum to a grid of 
NLTE synthetic spectra. 

Simultaneously fitting the U, B, V, Re, Ic, J, and H- 
bands light curves and radial velocity curve using the WD 
code, the geometrical and physical parameters of NSVS 1425 
were obtained. These results allow us to derive the absolute 
parameters of the system such as masses and radii of the 
components. 

We compare the position of the sdB star in NSVS 1425 
with other sdB and sdOB stars on the effective tempera- 
ture versus surface gravity diagram. We describe the possi- 
ble channels to form an sdB star in binaries and conclude 
that the post-common envelope development is probably the 
evolutionary scenario for NSVS 1425. The subsequent evo- 
lution of this system should lead to a cataclysmic variable. 
After a phase of angular momentum loss via gravitational 
radiation, this system will lie below the period gap of the 
cataclysmic variables. 
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